The finite mass beamsplitter in high power interferometers 
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The beamplitter in high-power interferometers is subject to significant radiation-pressure fluctuations. As 
a consequence, the phase relations which appear in the beamsplitter coupling equations oscillate and phase 
modulation fields are generated which add to the reflected fields. In this paper, the transfer function of the var- 
ious input fields impinging on the beamsplitter from all four ports onto the output field is presented including 
radiation-pressure effects. We apply the general solution of the coupling equations to evaluate the input-output 
relations of the dual-recycled laser-interferometer topology of the gravitational-wave detector GEO 600 and the 
power-recycling, signal-extraction topology of advanced LIGO. We show that the input-output relation exhibits 
a bright-port dark-port coupling. This mechanism is responsible for bright-port contributions to the noise den- 
sity of the output field and technical laser noise is expected to decrease the interferometer's sensitivity at low 
frequencies. It is shown quantitatively that the issue of technical laser noise is unimportant in this context if the 
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Earth-bound laser-interferometers seeking for gravitational waves dsn use high-power light fields in order to minimise 
the quantum-noise in the detection band. The main contribution to the quantum noise comes from the output port itself |5[. The 
output port vacuum field is reflected by the interferometer back towards the photodetector. Its noise spectral density was studied 
, in great detail One can manipulate the dark port field, thereby increasing the sensitivity of the gravitational-wave detection 

■ in the frequency band of interest. It was proposed by Caves to squeeze the vacuum field |8]. Combined with an appropriate 
ON ' filtering scheme,Jhe increase in sensitivity is limited by the squeezing factor. This was shown for different interferometer 
topologies |6, 9, 10]. In all these investigations radiation-pressure noise at the beamsplitter was not included. 

In this paper we analyze the effect of radiation-pressure fluctuations acting on the beamsplitter. We show that it gives rise 
to a coupling of the bright input port to the dark output port. Consequently this paper focusses on the field which enters the 
interferometer at the bright port, i. e. the input light which comes from the laser. We quantitatively investigate the contributions 
of quantum noise and technical laser noise of the bright input port to the interferometer's noise spectral density of the output 
field. We show that the coupling strongly depends on the interferometers topology and that technical laser noise might limit the 
detectors sensitivity at low frequencies. 

The paper is organized as follows. In section II, we discuss some general properties of beamsplitters and introduce the 
coupling equations of the fields. The solution of the coupling equations is presented in section III for a specific configuration, 
the power-recycled interferometer operating at dark fringe. In section IV, we present the noise spectral density for the current 
set-up of the dual-recycled GEO 600 interferometer and also for its envisioned design parameters. We conclude that in the design 
configuration of GEO 600 the dark port noise spectral density might be dominated by technical noise from the bright port at low 
frequencies. In section V, the same calculations are performed for the advanced LIGO configuration. A quantitative comparison 
shows that the relative contribution of the bright port noise to the output spectral density for advanced LIGO is smaller by five 
orders of magnitude than for GEO 600. 



II. THE COUPLING EQUATIONS 

Quantum fields are usually described by means of their annihilation and creation operators. The two-photon formalism devel- 
oped by Caves and Schumaker llHl turns out to be a more suitable formalism for measurements with heterodyne or homodyne 
detectors. These two classes of detectors measure the quadrature fields of the light whose amplitudes annihilate quanta of modu- 
lations. Correlations between the two sidebands built up by two-photon processes find a natural representation in that formalism 
and the spectral densities of the two quadratures' quantum noise is given by orthogonal sections through the so called noise 
ellipse. Modern publications discussing high-power interferometry show that one can derive simple and easy-to-interpret ex- 
pressions for the quadrature transfer functions of various configurations Therefore, we present all equations 
in the two-photon formalism benefitting from algebraic properties of the quadrature fields concerning radiation-pressure effects. 
The two quadrature amplitudes Si, a>2 merge into one single object which we call the quadrature vector: 
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FIG. 1: The noise density of the output field o determines the noise of the gravitational- wave detection. All c l and the input field i are 
propagating towards the beamsplitter. The fields d ! , o propagate away from it. The asymmetric beamsplitter reflects with a minus sign on the 
side where it is indicated in the picture. The A ! denote the classical amplitudes of the carrier in each direction. The south port does not contain 
any carrier field at dark fringe. 



Many important physical transformations acting on a can be interpreted geometrically as rotations and scalings of vectors in the 
space spanned by the two quadrature fields (i.e. the quadrature space). A more detailed treatment is given in 1 15]. Our notational 
conventions are introduced in Fig. ^ The classical carrier amplitudes are treated seperately from the modulation amplitudes. 
Therefore, we assume that the expectation values of both components of all quadrature vectors are much smaller than their 
carrier amplitudes A'. In the literature, one finds different conventions for the phase relations of the various fields which couple 
at the beamsplitter. However, for a lossless beamsplitter one can derive Stokes-like reciprocity relations involving the reflection 
and transmission of light which require that the amplitudes couple at the beamsplitter according to 

o = rc" - pP x c e 

d e = tc w - pP x i 

d" = t\ + P P- x c w 

d w = TC e +pP_ x c n (2) 

An explicit expression for P x is developed in the next section, when the power-recycled Michelson interferometer is discussed. 
At this point, we state some of its general features. The operator P x counts for the change of phase relations due to displacements 
x of the beamsplitter. Therefore, P x can be thought of as a propagator of the field along x corresponding to a rotation of the field's 
quadrature vector in quadrature space. Since the coupling equations relate modulation amplitudes, x denotes the amplitude for 
displacements of the beamsplitter at some frequency which is the modulation frequency of the field. Gravitational waves do not 
affect the motion of the beamsplitter in its own proper reference system and consequently x is independent of the gravitational- 
wave amplitude h. In that case, the equation of motion for x is completely determined by the radiation-pressure fluctuations of 
the light, i.e. by the fluctuations of the amplitude quadratures ai of all the fields. We derive the equation of motion in terms 
of momentum conservation. The beamsplitter has to compensate for the momentum flow of the ingoing and outgoing fields. 
Therefore, we make the following linearized ansatz in terms of the modulation amplitudes 

x oc A w (c^ + df ) - A e (cl + df) + A™(c™ + <). (3) 

The minus sign in front of the second bracket means that the momentum assigned to the east fields is carried in the opposite 
direction with respect to the momentum carried by the west fields, whereas the plus sign in front of the last bracket means that 
a motion of the beamsplitter downwards is equivalent to a motion towards the east concerning phase shifts of the reflected light. 
The square root of the spectral density of position fluctuations \J S{x) is supposed to be much smaller than the wavelength Ao 
of the carrier light. If that condition were not fulfilled, then the backaction of our measurement device on the test masses would 
be much larger than typical displacements induced by a gravitational wave (x « 10~ 10 • Ao, depending on the amplitude of 
the gravitational wave). If the field which the propagator P x acts on is not accompanied by a high-power carrier amplitude A, 
then the propagation becomes the unity matrix. In other words, the position fluctuations of the beamsplitter do not generate any 
sidebands, because there is no carrier on which sidebands with significant amplitude could be modulated. Then we obtain the 
following expression for the propagation along small displacements if the interferometer operates at dark fringe 

Pj = i 

P x c? = c> - A 3 ■ R(c\d l ) (4) 

The vector R is a linear function of x and consequently, it depends on all the fields which enter into the equation of motion of 
the beamsplitter. It should be clear that we just need one variable to determine the position of the beamsplitter since, concerning 
the phase shift of the reflected fields, a motion of the beamsplitter downwards is completely analogous to a motion to the right. 
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We are going to add three more equations to our system of coupling relations Eq. (0. The idea is to assign roundtrip 
transfer functions E, N, W and independent fields e, n, w to three of the four ports. The new fields comprise a sum of all fields 
originating in the corresponding port , e.g. vacuum fields due to losses or classical signal fields due to a gravitational wave. One 
may understand this step as some sort of closure of the ports by means of mirrors which reflect the outgoing light back to the 
beamsplitter. 

c e = Ed e +e 
-n = 7Vd" + fi 

c w = Wd w +w (5) 

The latter equations are the most general, linear equations which govern the roundtrip of the light. In the two-photon formalism, 
the transfer functions E, N, W are transfer matrices acting on quadrature vectors. 



III. THE INPUT-OUTPUT RELATION 



The input-output relations of an optical system comprise all contributions to the output field, i. e. the field which is detected 
by the photodiode. It is obtained by solving the coupling equations Eq. (|2) and Eq. (|5}: 

o = IO(i, w, n, e) (6) 

We present the solution of the coupling relations for a power-recycled interferometer with a 50/50 beamsplitter operating at dark 



FIG. 2: Power-recycled interferometer. Both arms are described by the same transfer matrix A. The west port contains a power-recycling 
mirror with amplitude reflectivity p pr . It forms the power-recycling cavity with the endmirrors of the Michelson interferometer. The mirror's 
distance to the beamsplitter is set to be an integer multiple of the carrier wavelength. The same holds for the pathlength of the light inside the 
two interferometer arms. 



fringe as shown in Fig. (0: 



1 50/50 beamsplitter 

A := A w = V2A e = V2A" 



P = T= -2 



A:=E = N 
W :=p pr (l + RPE) 



dark fringe condition 
power recycling 



(7) 



The power-recycling condition means that the transfer function W is a multiple of the identity map except for radiation-pressure 
effects. Without loss of generality, our condition requires that the distance of the beamsplitter to the power-recycling mirror is 
a multiple of the carrier wavelength, which also implies that the pathlength of the light inside the Michelson arms is a multiple 
of the carrier wavelength. In fact, the proper power-recycling condition is weaker than the one imposed here for simplicity. The 
proper condition merely requires that the combined pathlength through the Michelson arm and the power-recycling cavity is a 
multiple of the carrier wavelength. The radiation-pressure induced noise sidebands generated at the beamsplitter are derived 
from the matrix for small propagations and from the equation of motion of the beamsplitter. Small propagations P x lead to the 
following transformation of the quadrature vectors [15]: 



1 



^ A J / _ 
c UCcSd*) 



(8) 



Implicitly, we made use of the fact that c* is a modulation amplitude of a carrier field whose amplitude A J points in the direction 
of the amplitude quadrature of c? . The second term on the right-hand side corresponds to the noise sidebands which are excited 
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by fluctuations of the phase uiqx/c. In the two-photon formalism, phase fluctuations yield fluctuations of the phase quadrature 
whose noise amplitude is the phase shift multiplied by the amplitude of the carrier field. We assumed that the carrier frequency 
wo is much higher than the modulation frequency which, henceforth, is denoted by f2. The equation of motion for x is governed 
by Newton's law 



5P(c\d l ) 
mcfl 2 



(9) 



Here, Newton's equation is written in the domain of modulation frequencies and m is the mass of the beamsplitter. The fluctu- 
ating part SP of the radiation pressure is proportional to the right-hand side of Eq. (0. The constant of proportionality can be 
determined by comparing our expression for SP with expressions evaluated for simpler geometrical situations (e. g. see la]): 



SP = hw Q ■ A w (c^ + O - A e (cl + d\) + A n (c? + 



Bringing everything together, we cast Eq. (|8) into the form 
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(c-+d»)--i=(c e +d e ) + -i=(c" + d») 



(10) 



(11) 



where A denotes the amplitude of the light in the west port. Before we write down the input-ouput relation, we introduce the 
abbreviation 



K 




-K B 



K F 



fuv 2 A 2 
mc 2 Q 2 



(12) 



The coupling constant Kb is proportional to the power of the light at the beamsplitter by virtue of P = hw$A 2 . Inserting 
Eq. (II Q into the coupling equations and subsequently solving the system of linear equations for the output field, one obtains 



2(1 - A Ppi ) - (2 + (1 - A)p pi )(l + A)K 
(2A(1 - Ap w ) + (1 - (1 + 2p pr )A)(l + A)K)l 
+V2(1 - Ap pl )n 

-V2((l - Ap pi ) - (1 + p pi )(l + A)K)e 
+ (l + A) 2 Kw 



(13) 



If the radiation-pressure fluctuations are negligible, then the matrix K becomes zero and the input-output relations reduce to a 
well-known form. The most interesting aspect of this result is probably contained in the last term within the square brackets. 
It says that whenever there are radiation-pressure fluctuations acting on the beamsplitter, then fluctuations from the west port 
(also known as the bright port) couple to the output port. This contribution is proportional to the non-zero component of the 
matrix K. This might turn out to be a problem for all high-power interferometers, since the laser field suffers from high technical 
noise at low sideband frequencies, which couples into the field w. The technical noise at low frequencies can be several orders 
higher compared to pure vacuum fluctuations. The fact that the bright-port dark-port coupling is proportional to K also explains 
why the input-output relations are independent of radiation-pressure fluctuations acting on the power-recycling mirror. Those 
fluctuations are described by a matrix K 1 which has the same form than K and the transfer is governed by multiplying that 
matrix with K and K ■ K' is always zero. 



IV. THE NOISE SPECTRAL DENSITY OF THE GEO 600 TOPOLOGY 



In this section, we calculate the input-output relations of the dual-recycled configuration of GEO 600 and evaluate them in 
terms of the noise spectral density which is obtained under the following assumptions. The state of the input field i at the south 
port is a coherent vacuum field and w is the fraction of the laser field which transmits into the power-recycling cavity. Expressed 
in terms of single-sided spectral density matrices these properties assume the form 

S(i) = 1, S(W) = T 2 r ■ Stech(tt) (14) 

The matrix S tC ch is diagonal which means that our calculations do not account for correlations between the two quadratures built 
up inside the laser. The amount of technical noise which is brought into the interferometer by w is estimated from measurements 
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Signal-recycling 
mirror 



Detector Q 

FIG. 3: GEO 600 is a dual-recycled Michelson interferometer with a power- recycling mirror in the bright port that enhances the light power 
within the Michelson arms and a signal-recycling mirror in the dark port that can be tuned to a specific signal frequency. Since the arms are 
folded once, the effective armlengfh is doubled to 1200 m. The distance between the beampslitter and the so called far mirrors of the Michelson 
arms is 600 m, whereas the so called near mirrors which form the end of the arms are placed very close to the beamsplitter. 



performed on the GEO 600 laser. Optical losses occuring in real interferometers at the endmirrors or at the beamsplitter are not 
included in the sense that we do not mix the fields inside the interferometer with loss related vacuum fields. The value of the 
classical amplitude of the carrier light at different points of the interferometer is taken from real measurements. The equations 
of motion of all optical components are determined by the light pressure and the action of a gravitational wave. The latter one 
couples to the fields n and e. No significant signal is found in the field w since the distance of the power-recycling mirror to 
the origin of our reference frame (i. e the beamsplitter) is small compared to the lengths of the two Michelson arms. A transfer 
matrix A for the arms was first presented in 1 6] and was derived in |15] for the GEO 600 configuration applying the same 
formalism: 



A = ^{-K A l) ^ 

The optomechanical coupling constant K A of the Michelson arms is defined similarly to the beamsplitter coupling constant Kb 
in Eq. dl 2l i with the amplitude A substituted by the amplitude of the light inside the arms and the beamsplitter mass m substituted 
by the reduced mass for the two endmirrors (each having mass tum) which form the folded arms of GEO 600 (see Fig.|5Ji: 

^(A 2 /2) 

KA = 4 ' (m M /5)cW (16) 

By folding the arms, the effective armlength L becomes twice the distance between the far mirror and the beamsplitter. A 
gravitational wave h creates signal sidebands in both arms which possess equal amplitudes but different signs 

Ah 

(17) 



1 




(m M /5)Q 2 L 2 



The quantity h$QL is the standard quantum limit of GEO 600 with an infinite mass beamsplitter. The "true" quantum limit for 
GEO 600 also depends on the dynamics of the beamsplitter. We refrain from redefining h$QL in that manner, since here we want 
to discuss the beamsplitter dynamics explicitely and we do not want to find the reduced mass motion of the system. The problem 
to calculate the phase and coupling constant of a folded arm transfer function is related to the calculation of the same quantities 
for a delay line. A nice treatment of delay lines in our formalism can be found in the appendix of [ 13]. For this particular set of 
matrices [Eqs. (I15M17H . the input-output relation Eq. (1131 is given by 



i oy 2i 9± o o 

-Ki l) 1 + e \-K 2 



b + V2n (18) 



We substituted the field w by the transmitted bright port input field w = r pl b. The two constants K\ and K% depend on the arm 
and beamsplitter coupling constants 



Kt = K A + 2 cos' C-f)- K B , 

K 2 = 2cos 2 (^) -K B (19) 

1 - P pr e^ o 
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The coupling constant K2 is a product of 2 cos 2 (^) ~ 2 and the amplification factor for modulation fields inside the power- 
recycling cavity. We should emphasize that K 2 is independent of the arm coupling constant Ka and thus independent of the arm 
topology (i. e. whether it is a Michelson interferometer without arm cavities or with arm cavities). However, the modulus of the 
quantity K2 decreases if the arm length L is increased. From Eq. Jl 8I > one derives the input-output relation of the signal-recycled 
interferometer in the usual manner. Propagating fields from the beamsplitter to the signal-recycling mirror is accomplished by 
a rotation matrix D(<f>) acting in quadrature space which lacks the additional phase shift of the modulation fields since the 
wavelength A = (2ttc) /VL of the sidebands within the detection band (i. e. 10 Hz-1000 Hz) is much longer than the length of the 
signal-recycling cavity [7] 



yT/ V sm(0) cos(0) 



(20) 



The angle <fi is the detuning parameter of the signal-recycling cavity which is formed by the signal-recycling mirror and the 
Michelson interferometer. In Eq. Jl 81 . giving names Ti and Tf, to the transfer matrices of the fields i and b respectively, the 
input-output relation for the signal-recycled interferometer reads 



D{<t>)TiD{<!>) 
r D{0)-(T b b 



(p^TiDicfi - p sr 
fV2n 



1 i 



(21) 

The input-output relation determines the noise spectral density of the output field. The overall noise density is a sum of the 
two densities for the input field i sr and b. The latter one is the technical noise transferred from the bright port, the former 
one is the vacuum noise reflected at the dark port. It is convenient to normalize the spectral densities of the amplitude and 
phase quadratures of o sr such that the spectral density refers to the amplitude h of the gravitational wave which is contain in 
fi. The way how to do this normalization in matrix notation is shown in |9]. The evaluation of the spectral density is based 
on the parameter values according to Table [I] A detuning <j> = 0.015 means that the sideband which lies 600 Hz above the 





Symbol 


Value 


Light power at BS 


P 


300 W 


Transmissivity PRM 


r 2 


1.35% 


Transmissivity SRM 


r 2 

' sr 


2% 


Beamsplitter mass 


m 


9.3 kg 


Mirror mass 


m M 


5.6 kg 


Arm length 


L 


1200 m 


Frequency of laser 




1.77- 10 15 rad/s 


Detuning of SR cavity 


<t> 


0.015 rad 



TABLE I: Parameters of the GEO 600 configuration during the S3 run 1 18]. The detuning <f> of the signal-recycling cavity can be varied. The 
input light power at the power-recycling mirror was about 1.5W. 



carrier is resonantly amplified within the signal-recycling cavity. Adjusting the phase of the local oscillator in a homodyne 
detection scheme (corresponding to the electronic demodulation phase in heterodyne detection schemes), one can choose the 
direction in quadrature space along which the measurement is carried out. In that manner, the phase quadrature, the amplitude 
quadrature or some intermediate linear combination of these two can be measured. We refer to |16] for a deeper discussion 
of the quantum noise in heterodyne measurement schemes. Here, we restrict to measurements of the phase quadrature. The 
single-sided noise spectral density of the phase quadrature of the output field is shown in Fig.0] The bright port noise at low 
frequencies causes the optomechanical resonance to disappear from the noise spectral density. On the one hand, this effect is 
merely of theoretical interest as the currently measured noise density at low frequencies is dominated by seismic noise which 
couples to the optical fields through the mirror suspension. On the other hand, the result suggests that one has to investigate 
the role of bright port fluctuations for future interferometers. The beamsplitter coupling constant Kb is proportional to the light 
power at the beamsplitter. Therefore, one might expect that the transferred bright port noise becomes even more significant for 
high power interferometers of the next generation. The corresponding noise spectral density for GEO 600 with design power 
P = 10 kW and adjusted detuning tf> = 0.003 and transmissivity r s 2 r = 0.16% is shown in Fig. [5] assuming the same relative 
technical noise than before. Since the power of the carrier light is higher than in the previous case, all coupling constants are 
increased and the low frequencies noise experiences a shift upwards. Furthermore, the absolute technical bright-port noise was 
scaled by a factor 10 W/1.5 W derived from the two respective input powers. 

We conclude this section by suggesting a quantity which best characterizes the impact of the bright-port dark-port coupling 
on the output spectral density. That quantity should describe the balance of contributions coming from the bright port and the 
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Frequency [Hz] 



FIG. 4: Single-sided noise spectral density of GEO 600 with P = 300 W at the beamsplitter. The spectral density of the bright port vacuum 
field is lying below the dark port noise spectral density throughout the entire detection band. However, the technical noise from the bright 
port is dominating the spectral density up to 10 Hz where it is more than one order of magnitude higher than the vacuum noise density. The 
technical noise corresponds to an input laser field with power P m = 1.5 W . 




Frequency [Hz] 



FIG. 5: Single-sided spectral density of the dark-port field for the GEO 600 topology with design power P — 10 kW at the beamsplitter 
assuming the same relative technical noise here than for the low power interferometer underlying Fig.|4| However, the absolute technical noise 
of the input field is increased due to the higher input power of the light Pi n = 10 W. 



dark port to the output noise. We are looking for a characteristic function of the interferometer topology which is independent of 
the input power. We derive such a quantity from the input-output relation Eq. Jl 8i by comparing the components of the transfer 
matrix of the bright port field with the components of the matrix for the dark port field. At low frequencies (i. e. less than 10 Hz) 
it suffices to compare the values of the coupling constants K\ and K2 defined in Eq. dl9l . Their ratio \K\ \/\K^ | tells us which 
field mainly determines the fluctuations in the output field o. If the ratio is bigger than one, then the dark port field i dominates. 
If the ratio is less than one, then the bright port field b dominates. We call this ratio the low frequency balance and denote it by 
a. For GEO 600 without signal-recycling mirror one obtains the following expression: 

r pr / K A \ r pi rn + m M /5 

CTconv « -7T- 7777- + 1 = 7*— « 0-5 (22) 

2 \2K B J 2 m M /o 

This equation states that the bright-port fluctuations at low frequencies contribute twice as much as the dark-port fluctuations. 
Comparing the result with the noise spectral density of the vacuum fields in Fig.|4] we see that due to the signal-recycling mirror 
the bright port fluctuations become less important. The reason is that we approximate at frequencies which are less than the 
(half-)bandwidth of the signal-recycling cavity 7 sr = 200 Hz. For <C 2tt ■ 7 sr , the signal is weaker with signal-recycling 
cavity compared to a configuration without signal-recycling mirror. The bright port field behaves in the same way as the signal 
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field whereas the fluctuations from the dark port are nearly unaffected by the signal-recycling mirror. Therefore the noise-to- 
signal ratio of the dark-port fluctuations is increased with respect to the noise-to-signal ratio of the bright-port field which is not 
changed by the signal-recycling mirror. Since the detuning of the signal-recycling cavity is very small, we find the following 
balance for the dual-recycled configuration GEO 600 

_ r pr / K A \ r pr m + m M /5 

CTGEO » — 7777- + 1 = — 77— ~ 7-7 (23) 

r sr \2K B J r sr m M /o 

That value is in good agreement with the low frequency dark port and bright port vacuum noise in Fig.0] Inserting the respective 
values for the final set-up of GEO 600, the balance is <7fGEO ~ 15 which also agrees with the spectral densities in Fig. [5] 

V. THE NOISE SPECTRAL DENSITY FOR THE ADVANCED LIGO TOPOLOGY 

In this section, we apply the methods of the last section to the advanced LIGO configuration. According to the current plan 
IrHl . advanced LIGO will be a power-recycling, signal-extraction Michelson interferometer with arm cavities (see Fig.|6j. The 

i ETM 



PRM 



i ITM 



ITM ETM 



4 y'- h 



i SEM 

FIG. 6: The advanced LIGO configuration is a power-recycling, signal-extraction Michelson interferometer with arm cavities. The armlength 
is 4 km. Each arm cavity is formed by an input test mass and an end test mass. The signal-extraction cavity is formed by a mirror in the dark 
port and the interferometer which has the conventional LIGO topology. 

formulas which have to be applied for LIGO are identical to the formulas which were derived in the last section. The only 
difference lies in the definition of the arm coupling constant Ka and the signal field in terms of the standard quantum limit 
which counts for the arm topology of LIGO (compare with Eq. d!6t and Eq. tYTY ) 

8u; P 1 
A m M L 2 n 2 (j 2 + Sl 2 ) 



8h 
mML' 2 Q 2 



h SQL = — 7^77777 ( 25 ) 



Also the additional phase shift gained by the modulation fields which are now reflected at the inner test mass of the arm cavity 
has to be replaced by an expression which depends on the half-bandwidth 7 = cr 2 tra /AL of the arm cavity 

VLL /OA 

► arctan — (26) 



,7. 

All parameter values which enter the preceding definitions are gathered from 00 [3- They are listed in Tab.|H] The mass of 
the beamsplitter is accurate up to some small percentage. The noise spectral density of the output field shown in Fig.|7]lies well 
above the spectral density of the bright port vacuum and also above the technical noise from the bright port. Again the latter one 
is characterized by the same relative technical noise than in the two cases discussed for GEO 600. The fact that the bright-port 
dark-port coupling is insignificant for LIGO was anticipated and can be further quantified by performing a comparison of the 
beamsplitter and arm coupling constants. Due to the increased power in the arm cavities, LIGO's arm coupling constant Ka is 
much bigger than the beamsplitter coupling Kb. At low frequencies, the coupling constant Ka can be approximated by 

(4-)' (") 



mc 2 Q 2 \ T:? 
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Symbol 


Value 


Light power at BS 


P 


10 kW 


Transmissivity PRM 


r 2 


5% 


Transmissivity SRM 


T 2 
' sr 


19% 


Transmissivity of ITM 


2 

r itm 


3.3% 


Beamsplitter mass 


m 


13 kg 


Mirror mass 


m M 


40 kg 


Arm length 


L 


4000 m 


Frequency of laser 




1.77- 10 15 rad/s 


Detuning of SE cavity 


<t> 


tt/2 - 0.47 rad 



TABLE II: Parameters of the advanced LIGO configuration. Except for the beamsplitter mass, the values of the parameters are chosen 
according to Id I7 U1 711 . The transmissivity of the power-recycling mirror corresponds to a power amplification factor of 80 and so the input 
light power has to be 125 W. 




10 1 10 2 10 3 

Frequency [Hz] 



FIG. 7: Single-sided spectral density for the LIGO topology with P = 10 kW at the beamsplitter. We chose the same spectral density of 
relative technical laser noise than for the GEO 600 configurations. The bright-port fluctuations are negligible. 



Evaluating the ratio of the two matrix components Ki, Ki of Eq. (II 81 for a modulation frequency il — 2ir ■ 10 Hz which is far 
below the half -bandwidth of the arm cavities, we obtain 

2fl ( K A \ 8n m ( 4 \ 2 , 

tligo « hj#- + 1 » — ~ 8 • 10 4 - (28) 

r pr 7 \2K B J Tpr-yrriM \rf tm J 

Comparing with Eq. ( I23> . the result shows that the low frequency balance is five orders of magnitude bigger for LIGO than for 
GEO 600 corresponding to a weaker bright-port contribution to the output field. Even if the technical fluctuations of the input 
light are two or three orders of magnitude stronger than pure vacuum fluctuations, there will be no noticable contribution to the 
spectral density of the output field. 



VI. CONCLUSION 



We showed that the bright-port dark-port coupling gives rise to a significant contribution of technical fluctuations to the noise 
spectral density at low frequencies for a Michelson topology without arm cavities. In terms of the low frequency balance, we 
concluded that the LIGO topology exhibits a comparatively weak bright-port dark-port coupling relative to the contribution of 
the dark port noise. That is true even for a high level of technical laser noise. There are a couple of strategies to reduce these 
fluctuations at the dark port of the GEO 600 topology. One option is to decrease the transmissivity of the power-recycling mirror. 
The proposition seems to be in contradiction to Eq. d23t which states that the relative bright-port fluctuations increase with 
decreasing r pr . The reason why it works is, that the factor r pr /2 in front of the brackets has to be replaced by VLL/ (r pr c) if 
the following condition holds: r pr <C VtL/c. The required amplitude transmissivity had to be around lOppm which lies beyond 
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any practical feasibility and which is not desired for other reasons. The most obvious option is to increase the mass of the 
beamsplitter without increasing the masses of the endmirrors. One can see from Eq. fl23i that the beamsplitter mass has to be 
increased by two to three orders of magnitude depending on the technical noise which seems to be unfeasible again. The results 
of this paper do not have any practical relevance for currently operated interferometers whose sensitivity at low frequencies is 
limited not by the technical noise from the bright port, but either by seismic noise or thermal noise. The lasers which furnish 
the light for interferometers of the next generation are supposed to have a considerably lower amount of relative technical noise. 
Therefore, increasing the mass of the beamsplitter by a modest factor might already be sufficient in order to make the bright port 
fluctuations negligible for the GEO 600 topology. Also strategies which yield an increase of the sensitivity at high frequencies 
fioll are not influenced by our results. 

There is another generic mechanism existing by which a bright-port dark-port coupling is built up. If the transfer function of 
the two arms are not equal, then the input-output relation contains the following contribution from the bright port 

s = io ^ + i-£n + ~J%.w '* (29) 

where p and r denote the amplitude reflectivity and transmissivity of the beamsplitter. There are different reasons why the two 
arms are not described by the same transfer function. One reason could be that one arm is detuned intentionally in order to 
transmit some carrier light towards the detector where it serves as local oscillator for a homodyne measurement of the signal. 
For a small detuning of just one arm, the carrier light is transmitted into the phase quadrature of o governed by a transfer function 
which is proportional to the detuning. If the carrier light is to be transmitted into the amplitude quadrature of 5, then one can 
choose the detuning of one arm to be tt/2 and the detuning of the other arm to be slightly less than tt/2. In that case, the west 
port has to be detuned accordingly to maintain the power-recycling condition (which cannot hold exactly here, since the two 
arms are detuned relative to each other). An unintentional reason for different transfer functions of the two arms could be that 
the transmissivity and the reflectivity of the beamsplitter are not the same. Then, radiation-pressure fluctuations in the two arms 
would be different by virtue of the different powers of the two respective carrier fields. Carrier light is then transmitted into 
the phase quadrature of o and the corresponding transfer function is proportional to the difference of the power reflectivity and 
transmissivity of the beamsplitter. There would also be a small loss of the optical signal due to a partial transmission into the 
bright port which is proportional to the same difference of power reflectivity and transmissivity. 
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The beamplitter in high-power interferometers is subject to significant radiation-pressure fluctuations. As 
a consequence, the phase relations which appear in the beamsplitter coupling equations oscillate and phase 
modulation fields are generated which add to the reflected fields. In this paper, the transfer function of the var- 
ious input fields impinging on the beamsplitter from all four ports onto the output field is presented including 
radiation-pressure effects. We apply the general solution of the coupling equations to evaluate the input-output 
relations of the dual-recycled laser-interferometer topology of the gravitational-wave detector GEO 600 and the 
power-recycling, signal-extraction topology of advanced LIGO. We show that the input-output relation exhibits 
a bright-port dark-port coupling. This mechanism is responsible for bright-port contributions to the noise den- 
sity of the output field and technical laser noise is expected to decrease the interferometer's sensitivity at low 
frequencies. It is shown quantitatively that the issue of technical laser noise is unimportant in this context if the 
interferometer contains arm cavities. 



PACS numbers: 04.80.Nn, 03.65.Ta, 42.50.Dv, 95.55.Ym 

I. INTRODUCTION 

Earth-bound laser-interferometers seeking for gravitational 
waves [1-4] use high-power light fields in order to minimise 
the quantum-noise in the detection band. The main contribu- 
tion to the quantum noise comes from the output port itself [5] . 
The output port vacuum field is reflected by the interferometer 
back towards the photodetector. Its noise spectral density was 
studied in great detail [6, 7] . One can manipulate the dark port 
field, thereby increasing the sensitivity of the gravitational- 
wave detection in the frequency band of interest. It was pro- 
posed by Caves to squeeze the vacuum field [8]. Combined 
with an appropriate filtering scheme, the increase in sensitiv- 
ity is limited by the squeezing factor. This was shown for 
different interferometer topologies [6, 9, 10]. In all these in- 
vestigations radiation-pressure noise at the beamsplitter was 
not included. 

In this paper we analyze the effect of radiation-pressure 
fluctuations acting on the beamsplitter. We show that it gives 
rise to a coupling of the bright input port to the dark output 
port. Consequently this paper focusses on the field which en- 
ters the interferometer at the bright port, i. e. the input light 
which comes from the laser. We quantitatively investigate the 
contributions of quantum noise and technical laser noise of the 
bright input port to the interferometer's noise spectral density 
of the output field. We show that the coupling strongly de- 
pends on the interferometers topology and that technical laser 
noise might limit the detectors sensitivity at low frequencies. 

The paper is organized as follows. In section II, we dis- 
cuss some general properties of beamsplitters and introduce 
the coupling equations of the fields. The solution of the cou- 
pling equations is presented in section III for a specific config- 
uration, the power-recycled interferometer operating at dark 
fringe. In section IV, we present the noise spectral density for 
the current set-up of the dual-recycled GEO 600 interferom- 
eter and also for its envisioned design parameters. We con- 
clude that in the design configuration of GEO 600 the dark 
port noise spectral density might be dominated by technical 



noise from the bright port at low frequencies. In section V, the 
same calculations are performed for the advanced LIGO con- 
figuration. A quantitative comparison shows that the relative 
contribution of the bright port noise to the output spectral den- 
sity for advanced LIGO is smaller by five orders of magnitude 
than for GEO 600. 



II. THE COUPLING EQUATIONS 

Quantum fields are usually described by means of their an- 
nihilation and creation operators. The two-photon formalism 
developed by Caves and Schumaker [11] turns out to be a 
more suitable formalism for measurements with heterodyne 
or homodyne detectors. These two classes of detectors mea- 
sure the quadrature fields of the light whose amplitudes anni- 
hilate quanta of modulations. Correlations between the two 
sidebands built up by two-photon processes find a natural rep- 
resentation in that formalism and the spectral densities of the 
two quadratures' quantum noise is given by orthogonal sec- 
tions through the so called noise ellipse. Modern publications 
discussing high-power interferometry show that one can de- 
rive simple and easy-to-interpret expressions for the quadra- 
ture transfer functions of various configurations [6, 7, 9, 12- 
14]. Therefore, we present all equations in the two-photon for- 
malism benefitting from algebraic properties of the quadrature 
fields concerning radiation-pressure effects. The two quadra- 
ture amplitudes &i, a 2 merge into one single object which we 
call the quadrature vector: 




Many important physical transformations acting on a can be 
interpreted geometrically as rotations and scalings of vectors 
in the space spanned by the two quadrature fields (i.e. the 
quadrature space). A more detailed treatment is given in [15]. 
Our notational conventions are introduced in Fig. 1 . The clas- 
sical carrier amplitudes are treated seperately from the modu- 
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FIG. 1: The noise density of the output field o determines the noise 
of the gravitational-wave detection. All c l and themput field i are 
propagating towards the beamsplitter. The fields d\ o propagate 
away from it. The asymmetric beamsplitter reflects with a minus 
sign on the side where it is indicated in the picture. The A 1 denote 
the classical amplitudes of the carrier in each direction. The south 
port does not contain any carrier field at dark fringe. 

lation amplitudes. Therefore, we assume that the expectation 
values of both components of all quadrature vectors are much 
smaller than their carrier amplitudes A 1 . In the literature, one 
finds different conventions for the phase relations of the var- 
ious fields which couple at the beamsplitter. However, for a 
lossless beamsplitter one can derive Stokes-like reciprocity re- 
lations involving the reflection and transmission of light which 
require that the amplitudes couple at the beamsplitter accord- 
ing to 

O = TC" - pP x C e 

d e = tc w - pP~i 

d" = t\ + p P_ x c w 

d w = rc e +pP- x c n (2) 

An explicit expression for P x is developed in the next sec- 
tion, when the power-recycled Michelson interferometer is 
discussed. At this point, we state some of its general features. 
The operator P x counts for the change of phase relations due 
to displacements x of the beamsplitter. Therefore, P x can be 
thought of as a propagator of the field along x correspond- 
ing to a rotation of the field's quadrature vector in quadra- 
ture space. Since the coupling equations relate modulation 
amplitudes, x denotes the amplitude for displacements of the 
beamsplitter at some frequency f2 which is the modulation fre- 
quency of the field. Gravitational waves do not affect the mo- 
tion of the beamsplitter in its own proper reference system 
and consequently x is independent of the gravitational-wave 
amplitude h. In that case, the equation of motion for x, is 
completely determined by the radiation-pressure fluctuations 
of the light, i.e. by the fluctuations of the amplitude quadra- 
tures a\ of all the fields. We derive the equation of motion 
in terms of momentum conservation. The beamsplitter has to 
compensate for the momentum flow of the ingoing and outgo- 
ing fields. Therefore, we make the following linearized ansatz 
in terms of the modulation amplitudes 

x oc A w (c? + <) - A e (c^ + d\) + A"(cl l + <). (3) 

The minus sign in front of the second bracket means that the 
momentum assigned to the east fields is carried in the oppo- 



site direction with respect to the momentum carried by the 
west fields, whereas the plus sign in front of the last bracket 
means that a motion of the beamsplitter downwards is equiv- 
alent to a motion towards the east concerning phase shifts of 
the reflected light. The square root of the spectral density of 
position fluctuations \J S(x) is supposed to be much smaller 
than the wavelength Ao of the carrier light. If that condition 
were not fulfilled, then the backaction of our measurement de- 
vice on the test masses would be much larger than typical dis- 
placements induced by a gravitational wave (x w 10~ 10 • Ao, 
depending on the amplitude of the gravitational wave). If the 
field which the propagator P x acts on is not accompanied by 
a high-power carrier amplitude A, then the propagation be- 
comes the unity matrix. In other words, the position fluctu- 
ations of the beamsplitter do not generate any sidebands, be- 
cause there is no carrier on which sidebands with significant 
amplitude could be modulated. Then we obtain the following 
expression for the propagation along small displacements if 
the interferometer operates at dark fringe 

Pj = i 

P x c? = c j - A 3 • S(c\d*) (4) 

The vector R is a linear function of x and consequently, it de- 
pends on all the fields which enter into the equation of motion 
of the beamsplitter. It should be clear that we just need one 
variable to determine the position of the beamsplitter since, 
concerning the phase shift of the reflected fields, a motion of 
the beamsplitter downwards is completely analogous to a mo- 
tion to the right. 

We are going to add three more equations to our system 
of coupling relations Eq. (2). The idea is to assign roundtrip 
transfer functions E, N, W and independent fields e, n, w to 
three of the four ports. The new fields comprise a sum of 
all fields originating in the corresponding port , e.g. vacuum 
fields due to losses or classical signal fields due to a gravita- 
tional wave. One may understand this step as some sort of 
closure of the ports by means of mirrors which reflect the out- 
going light back to the beamsplitter. 

c e = Ed e + e 
c" = Afd n + n 

c w = Wd w +w (5) 

The latter equations are the most general, linear equations 
which govern the roundtrip of the light. In the two-photon 
formalism, the transfer functions E, N, W are transfer matri- 
ces acting on quadrature vectors. 



III. THE INPUT-OUTPUT RELATION 

The input-output relations of an optical system comprise 
all contributions to the output field, i. e. the field which is 
detected by the photodiode. It is obtained by solving the cou- 
pling equations Eq. (2) and Eq. (5): 

o = IO(i, w, n, e) (6) 
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FIG. 2: Power-recycled interferometer. Both arms are described by 
the same transfer matrix A. The west port contains a power-recycling 
mirror with amplitude reflectivity p pr . It forms the power-recycling 
cavity with the endmirrors of the Michelson interferometer. The mir- 
ror's distance to the beamsplitter is set to be an integer multiple of 
the carrier wavelength. The same holds for the pathlength of the 
light inside the two interferometer arms. 



We present the solution of the coupling relations for a power- 
recycled interferometer with a 50/50 beamsplitter operating at 
dark fringe as shown in Fig. (2): 



P = T =—2 

A:=A W = V2A e = V2A™ 
A := E = N 
W:=p pr (l + RPE) 



50/50 beamsplitter 

dark fringe condition 
power recycling (7) 



The power-recycling condition means that the transfer func- 
tion W is a multiple of the identity map except for radiation- 
pressure effects. Without loss of generality, our condition 
requires that the distance of the beamsplitter to the power- 
recycling mirror is a multiple of the carrier wavelength, which 
also implies that the pathlength of the light inside the Michel- 
son arms is a multiple of the carrier wavelength. In fact, the 
proper power-recycling condition is weaker than the one im- 
posed here for simplicity. The proper condition merely re- 
quires that the combined pathlength through the Michelson 
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arm and the power-recycling cavity is a multiple of the car- 
rier wavelength. The radiation-pressure induced noise side- 
bands generated at the beamsplitter are derived from the ma- 
trix for small propagations and from the equation of motion of 
the beamsplitter. Small propagations P x lead to the following 
transformation of the quadrature vectors [15]: 
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^(cScr) 



(8) 



Implicitly, we made use of the fact that c J is a modulation am- 
plitude of a carrier field whose amplitude A 3 points in the di- 
rection of the amplitude quadrature of cP . The second term on 
the right-hand side corresponds to the noise sidebands which 
are excited by fluctuations of the phase ojqx/c. In the two- 
photon formalism, phase fluctuations yield fluctuations of the 
phase quadrature whose noise amplitude is the phase shift 
multiplied by the amplitude of the carrier field. We assumed 
that the carrier frequency uj is much higher than the mod- 
ulation frequency which, henceforth, is denoted by fL The 
equation of motion for x is governed by Newton's law 



x — — - 



5P{c\d l ) 
mcfl 2 



(9) 



Here, Newton's equation is written in the domain of modula- 
tion frequencies and m is the mass of the beamsplitter. The 
fluctuating part SP of the radiation pressure is proportional 
to the right-hand side of Eq. (3). The constant of proportion- 
ality can be determined by comparing our expression for 5P 
with expressions evaluated for simpler geometrical situations 
(e. g. see [6]): 



SP = hw ■ [A- W {cf + O - K e {c\ + d\) + A n (c? + d?) 

(10) 

Bringing everything together, we cast Eq. (8) into the form 



hulAA j (0 



mc 



(c»+d-)--L(c e +d e ) + -L(c" + d«) 



(11) 
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where A denotes the amplitude of the light in the west port. 
Before we write down the input-ouput relation, we introduce 
the abbreviation 



the system of linear equations for the output field, one obtains 



K = 





-Kb 



K, 



hu 2 A 2 



mc 



2 n 2 



(12) 



The coupling constant K B is proportional to the power of the 
light at the beamsplitter by virtue of P = hw h 2 . Inserting 
Eq. (11) into the coupling equations and subsequently solving 



2(1 - Ap pr ) - (2 + (1 - A) PpT ){l + A)K 
{2A{1 - A Pw ) + (1 - (1 + 2p pi )A)(l + A)K)\ 
+V2(1 - A Pw )n 

-V2((l - A Ppr ) - (1 + p pr )(l + A)K)e 



+ (1 + A) 2 Kw 



(13) 
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If the radiation-pressure fluctuations are negligible, then the 
matrix K becomes zero and the input-output relations reduce 
to a well-known form. The most interesting aspect of this re- 
sult is probably contained in the last term within the square 
brackets. It says that whenever there are radiation-pressure 
fluctuations acting on the beamsplitter, then fluctuations from 
the west port (also known as the bright port) couple to the 
output port. This contribution is proportional to the non-zero 
component of the matrix K. This might turn out to be a prob- 
lem for all high-power interferometers, since the laser field 
suffers from high technical noise at low sideband frequencies, 
which couples into the field w. The technical noise at low 
frequencies can be several orders higher compared to pure 
vacuum fluctuations. The fact that the bright-port dark-port 
coupling is proportional to K also explains why the input- 
output relations are independent of radiation-pressure fluctua- 
tions acting on the power-recycling mirror. Those fluctuations 
are described by a matrix K' which has the same form than K 
and the transfer is governed by multiplying that matrix with 
K and K ■ K' is always zero. 



IV. THE NOISE SPECTRAL DENSITY OF THE GEO 600 
TOPOLOGY 

In this section, we calculate the input-output relations of the 
dual-recycled configuration of GEO 600 and evaluate them in 
terms of the noise spectral density which is obtained_ under 
the following assumptions. The state of the input field i at the 
south port is a coherent vacuum field and w is the fraction of 
the laser field which transmits into the power-recycling cavity. 
Expressed in terms of single-sided spectral density matrices 
these properties assume the form 



S(i) = l, S(w)=r, ; 



'tech 



(fi) 



(14) 



The matrix Stcch is diagonal which means that our calcula- 
tions do not account for correlations between the two quadra- 
tures built up inside the laser. The amount of technical noise 
which is brought into the interferometer by w is estimated 
from measurements performed on the GEO 600 laser. Optical 
losses occuring in real interferometers at the endmirrors or at 
the beamsplitter are not included in the sense that we do not 
mix the fields inside the interferometer with loss related vac- 
uum fields. The value of the classical amplitude of the carrier 
light at different points of the interferometer is taken from real 
measurements. The equations of motion of all optical compo- 
nents are determined by the light pressure and the action of 
a gravitational wave. The latter one couples to the fields n 
and e. No significant signal is found in the field w since the 
distance of the power-recycling mirror to the origin of our ref- 
erence frame (i. e the beamsplitter) is small compared to the 
lengths of the two Michelson arms. A transfer matrix A for 
the arms was first presented in [6] and was derived in [15] for 
the GEO 600 configuration applying the same formalism: 



1 
-K A 1 



(15) 




Detector Q 



FIG. 3: GEO 600 is a dual-recycled Michelson interferometer with 
a power-recycling mirror in the bright port that enhances the light 
power within the Michelson arms and a signal-recycling mirror in the 
dark port that can be tuned to a specific signal frequency. Since the 
arms are folded once, the effective armlength is doubled to 1200 m. 
The distance between the beampslitter and the so called far mirrors 
of the Michelson arms is 600 m, whereas the so called near mirrors 
which form the end of the arms are placed very close to the beam- 
splitter. 



The optomechanical coupling constant Ka of the Michelson 
arms is defined similarly to the beamsplitter coupling constant 
Kb in Eq. (12) with the amplitude A substituted by the ampli- 
tude of the light inside the arms and the beamsplitter mass m 
substituted by the reduced mass for the two endmirrors (each 
having mass Tom) which form the folded arms of GEO 600 
(see Fig. 3): 



K A = A- 



/^g(A 2 /2) 
(m M /5)c 2 ft 2 



(16) 



By folding the arms, the effective armlength L becomes twice 
the distance between the far mirror and the beamsplitter. A 
gravitational wave h creates signal sidebands in both arms 
which possess equal amplitudes but different signs 



n 



^-SQL 



h 2 

"SQL 



5H 



m M 



n 2 L 2 



(17) 



The quantity Hsql is the standard quantum limit of GEO 600 
with an infinite mass beamsplitter. The "true" quantum limit 
for GEO 600 also depends on the dynamics of the beamsplit- 
ter. We refrain from redefining hsQL in that manner, since 
here we want to discuss the beamsplitter dynamics explicitely 
and we do not want to find the reduced mass motion of the 
system. The problem to calculate the phase and coupling con- 
stant of a folded arm transfer function is related to the calcu- 
lation of the same quantities for a delay line. A nice treatment 
of delay lines in our formalism can be found in the appendix 
of [13]. For this particular set of matrices [Eqs. (15), (17)], the 
input-output relation Eq. (13) is given by 



n 1i — 



1 

-K x 1 



OL f 

-K 2 



b + V2ri (18) 
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We substituted the field w by the transmitted bright port input 
field w = T pr b. The two constants K\ and K 2 depend on the 
arm and beamsplitter coupling constants 



K x = K A + 2cos 2 ( l -f)-K B , 



K 2 



2cos a (^) 



' pr 



1 - p pr e 2 



(19) 



The coupling constant K 2 is a product of 2cos 2 (^) w 2 
and the amplification factor for modulation fields inside the 
power-recycling cavity. We should emphasize that K 2 is in- 
dependent of the arm coupling constant Ka and thus inde- 
pendent of the arm topology (i. e. whether it is a Michel- 
son interferometer without arm cavities or with arm cavities). 
However, the modulus of the quantity K 2 decreases if the 
arm length L is increased. From Eq. (18) one derives the 
input-output relation of the signal-recycled interferometer in 
the usual manner. Propagating fields from the beamsplitter 
to the signal-recycling mirror is accomplished by a rotation 
matrix D(<f>) acting in quadrature space which lacks the ad- 
ditional phase shift of the modulation fields since the wave- 
length A = (27rc)/0 of the sidebands within the detection 
band (i. e. 10 Hz-1000 Hz) is much longer than the length of 
the signal-recycling cavity [7] 



u(0 = (W) ~ si ^) 

v ' Vsm(0) cos(0) 



(20) 



The angle (j) is the detuning parameter of the signal-recycling 
cavity which is formed by the signal-recycling mirror and the 
Michelson interferometer. In Eq. (18), giving names Tj and 
Th to the transfer matrices of the fields i and b respectively, 
the input-output relation for the signal-recycled interferometer 
reads 



\-p sr - D{4,)TiD{cj)) 

+T sr D(4>) ■ (T b b + V2nj 



(pWTiDtf) - Ps 



(21) 



The input-output relation determines the noise spectral den- 
sity of the output field. The overall noise density is a sum of 
the two densities for the input field i sr and b. The latter one is 
the technical noise transferred from the bright port, the former 
one is the vacuum noise reflected at the dark port. It is con- 
venient to normalize the spectral densities of the amplitude 
and phase quadratures of o sr such that the spectral density 
refers to the amplitude h of the gravitational wave which is 
contain in n. The way how to do this normalization in ma- 
trix notation is shown in [9]. The evaluation of the spectral 
density is based on the parameter values according to Table 
I. A detuning <j> = 0.015 means that the sideband which 
lies 600 Hz above the carrier is resonantly amplified within 
the signal-recycling cavity. Adjusting the phase of the lo- 
cal oscillator in a homodyne detection scheme (corresponding 
to the electronic demodulation phase in heterodyne detection 
schemes), one can choose the direction in quadrature space 
along which the measurement is carried out. In that manner, 
the phase quadrature, the amplitude quadrature or some inter- 
mediate linear combination of these two can be measured. We 





Symbol 


Value 


Light power at BS 


P 


300 W 


Transmissivity PRM 


r 2 

1 pr 


1.35% 


Transmissivity SRM 




2% 


Beamsplitter mass 


m 


9.3 kg 


Mirror mass 


m M 


5.6 kg 


Arm length 


L 


1200 m 


Frequency of laser 




1.77- 10 15 rad/s 


Detuning of SR cavity 


<f> 


0.015 rad 



TABLE I: Parameters of the GEO 600 configuration during the S3 
run [18]. The detuning <f> of the signal -recycling cavity can be varied. 
The input light power at the power-recycling mirror was about 1.5W. 



refer to [16] for a deeper discussion of the quantum noise in 
heterodyne measurement schemes. Here, we restrict to mea- 
surements of the phase quadrature. The single-sided noise 
spectral density of the phase quadrature of the output field 
is shown in Fig. 4. The bright port noise at low frequencies 
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FIG. 4: Single-sided noise spectral density of GEO 600 with P = 
300 W at the beamsplitter. The spectral density of the bright port vac- 
uum field is lying below the dark port noise spectral density through- 
out the entire detection band. However, the technical noise from the 
bright port is dominating the spectral density up to 10 Hz where it 
is more than one order of magnitude higher than the vacuum noise 
density. The technical noise corresponds to an input laser field with 
power Pi n = 1.5 W . 

causes the optomechanical resonance to disappear from the 
noise spectral density. On the one hand, this effect is merely 
of theoretical interest as the currently measured noise density 
at low frequencies is dominated by seismic noise which cou- 
ples to the optical fields through the mirror suspension. On 
the other hand, the result suggests that one has to investigate 
the role of bright port fluctuations for future interferometers. 
The beamsplitter coupling constant Kb is proportional to the 
light power at the beamsplitter. Therefore, one might expect 
that the transferred bright port noise becomes even more sig- 
nificant for high power interferometers of the next generation. 



6 



The corresponding noise spectral density for GEO 600 with 
design power P = 10 kW and adjusted detuning <f> = 0.003 
and transmissivity r s 2 r = 0.16% is shown in Fig. 5 assum- 
ing the same relative technical noise than before. Since the 




Frequency [Hz] 

FIG. 5: Single-sided spectral density of the dark-port field for the 
GEO 600 topology with design power P = 10 kW at the beam- 
splitter assuming the same relative technical noise here than for the 
low power interferometer underlying Fig. 4. However, the absolute 
technical noise of the input field is increased due to the higher input 
power of the light P in = 10 W. 



tuations. Comparing the result with the noise spectral den- 
sity of the vacuum fields in Fig. 4, we see that due to the 
signal-recycling mirror the bright port fluctuations become 
less important. The reason is that we approximate at frequen- 
cies which are less than the (half-)bandwidth of the signal- 
recycling cavity 7 sr = 200 Hz. For <C 2tt ■ j sr , the signal 
is weaker with signal-recycling cavity compared to a config- 
uration without signal-recycling mirror. The bright port field 
behaves in the same way as the signal field whereas the fluctu- 
ations from the dark port are nearly unaffected by the signal- 
recycling mirror. Therefore the noise-to-signal ratio of the 
dark-port fluctuations is increased with respect to the noise- 
to- signal ratio of the bright-port field which is not changed by 
the signal-recycling mirror. Since the detuning of the signal- 
recycling cavity is very small, we find the following balance 
for the dual-recycled configuration GEO 600 

T" pr / ^ \ r pr to + tom/5 
CTGEO ~ 7777- + 1 = 77- ~ 7 ■< ( 23 ) 

r sr \2K B J r sr m M / 5 

That value is in good agreement with the low frequency dark 
port and bright port vacuum noise in Fig. 4. Inserting the re- 
spective values for the final set-up of GEO 600, the balance is 
ffGEO ~ 15 which also agrees with the spectral densities in 
Fig. 5. 

V. THE NOISE SPECTRAL DENSITY FOR THE 
ADVANCED LIGO TOPOLOGY 



power of the carrier light is higher than in the previous case, 
all coupling constants are increased and the low frequencies 
noise experiences a shift upwards. Furthermore, the absolute 
technical bright-port noise was scaled by a factor 10 W/1.5 W 
derived from the two respective input powers. 

We conclude this section by suggesting a quantity which 
best characterizes the impact of the bright-port dark-port cou- 
pling on the output spectral density. That quantity should 
describe the balance of contributions coming from the bright 
port and the dark port to the output noise. We are looking for 
a characteristic function of the interferometer topology which 
is independent of the input power. We derive such a quantity 
from the input-output relation Eq. (18) by comparing the com- 
ponents of the transfer matrix of the bright port field with the 
components of the matrix for the dark port field. At low fre- 
quencies (i. e. less than 10 Hz) it suffices to compare the val- 
ues of the coupling constants K\ and K2 defined in Eq. (19). 
Their ratio | K\ | / 1 K2 | tells us which field mainly determines 
the fluctuations in the output field o. If the ratio is bigger 
than one, then the dark port field i dominates. If the ratio is 
less than one, then the bright port field b dominates. We call 
this ratio the low frequency balance and denote it by a. For 
GEO 600 without signal-recycling mirror one obtains the fol- 
lowing expression: 



/ K A 
2 \2K B 



1 pr 



1 ] = ^ 



r pr to + tom/5 



to m / 5 



0.5 (22) 



This equation states that the bright-port fluctuations at low 
frequencies contribute twice as much as the dark-port fluc- 



In this section, we apply the methods of the last section to 
the advanced LIGO configuration. According to the current 
plan [17], advanced LIGO will be a power-recycling, signal- 
extraction Michelson interferometer with arm cavities (see 
Fig. 6). The formulas which have to be applied for LIGO are 



' ETM 



PRM 



' ITM 
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ETM 



H 7'- h 



1 SEM 



FIG. 6: The advanced LIGO configuration is a power-recycling, 
signal-extraction Michelson interferometer with arm cavities. The 
armlength is 4 km. Each arm cavity is formed by an input test mass 
and an end test mass. The signal-extraction cavity is formed by a 
mirror in the dark port and the interferometer which has the conven- 
tional LIGO topology. 

identical to the formulas which were derived in the last sec- 
tion. The only difference lies in the definition of the arm cou- 
pling constant Ka and the signal field in terms of the standard 
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quantum limit which counts for the arm topology of LIGO 
(compare with Eq. (16) and Eq. (17)) 



10" 



K A 

h-SQL 



8uj a P 



1 



m M L 2 ft 2 ( 7 2 + ft 2 ) 
8h 



(24) 
(25) 



Also the additional phase shift gained by the modulation fields 
which are now reflected at the inner test mass of the arm cavity 
has to be replaced by an expression which depends on the half- 
bandwidth 7 = cr 2 tm /4L of the arm cavity 



m /ft 

> arctan — 

c \7 



(26) 



All parameter values which enter the preceding definitions are 
gathered from [6, 7, 17]. They are listed in Tab. II. The mass 
of the beamsplitter is accurate up to some small percentage. 
The noise spectral density of the output field shown in Fig. 7 





Symbol 




Value 


Light power at BS 


P 




10 kW 


Transmissivity PRM 


2 




5% 


Transmissivity SRM 


r 2 

' sr 




19% 


Transmissivity of ITM 






3.3% 


Beamsplitter mass 


m 




13 kg 


Mirror mass 


m M 




40 kg 


Arm length 


L 




4000 m 


Frequency of laser 


wo 


1.77- 


10 15 rad/s 


Detuning of SE cavity 


<t> 


tt/2 


- 0.47 rad 



TABLE II: Parameters of the advanced LIGO configuration. Except 
for the beamsplitter mass, the values of the parameters are chosen 
according to [6, 7, 17]. The transmissivity of the power-recycling 
mirror corresponds to a power amplification factor of 80 and so the 
input light power has to be 125 W. 

lies well above the spectral density of the bright port vacuum 
and also above the technical noise from the bright port. Again 
the latter one is characterized by the same relative techni- 
cal noise than in the two cases discussed for GEO 600. The 
fact that the bright-port dark-port coupling is insignificant for 
LIGO was anticipated and can be further quantified by per- 
forming a comparison of the beamsplitter and arm coupling 
constants. Due to the increased power in the arm cavities, 
LIGO's arm coupling constant Ka is much bigger than the 
beamsplitter coupling Kb- At low frequencies, the coupling 
constant K A can be approximated by 



K A 



8uj P 
HP 



me- 



al) 



Evaluating the ratio of the two matrix components K\ , K2 of 
Eq. (18) for a modulation frequency ft = 2-7T • 10 Hz which is 
far below the half -bandwidth of the arm cavities, we obtain 



cligo 



2ft / K A 



Tp r 7 \2K B 



1 



8ft 



Tpr^niM \T itm 



(28) 
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FIG. 7: Single-sided spectral density for the LIGO topology with 
P — 10 kW at the beamsplitter. We chose the same spectral density 
of relative technical laser noise than for the GEO 600 configurations. 
The bright-port fluctuations are negligible. 



Comparing with Eq. (23), the result shows that the low fre- 
quency balance is five orders of magnitude bigger for LIGO 
than for GEO 600 corresponding to a weaker bright-port con- 
tribution to the output field. Even if the technical fluctuations 
of the input light are two or three orders of magnitude stronger 
than pure vacuum fluctuations, there will be no noticable con- 
tribution to the spectral density of the output field. 



VI. CONCLUSION 

We showed that the bright-port dark-port coupling gives 
rise to a significant contribution of technical fluctuations to 
the noise spectral density at low frequencies for a Michel- 
son topology without arm cavities. In terms of the low fre- 
quency balance, we concluded that the LIGO topology ex- 
hibits a comparatively weak bright-port dark-port coupling 
relative to the contribution of the dark port noise. That is 
true even for a high level of technical laser noise. There are 
a couple of strategies to reduce these fluctuations at the dark 
port of the GEO 600 topology. One option is to decrease the 
transmissivity of the power-recycling mirror. The proposition 
seems to be in contradiction to Eq. (23) which states that the 
relative bright-port fluctuations increase with decreasing r pr . 
The reason why it works is, that the factor r pr /2 in front of 
the brackets has to be replaced by ftL/(r pr c) if the follow- 
ing condition holds: Tp r <C ft-L/c. The required amplitude 
transmissivity had to be around lOppm which lies beyond any 
practical feasibility and which is not desired for other reasons. 
The most obvious option is to increase the mass of the beam- 
splitter without increasing the masses of the endmirrors. One 
can see from Eq. (23) that the beamsplitter mass has to be in- 
creased by two to three orders of magnitude depending on the 
technical noise which seems to be unfeasible again. The re- 
sults of this paper do not have any practical relevance for cur- 
rently operated interferometers whose sensitivity at low fre- 
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quencies is limited not by the technical noise from the bright 
port, but either by seismic noise or thermal noise. The lasers 
which furnish the light for interferometers of the next gener- 
ation are supposed to have a considerably lower amount of 
relative technical noise. Therefore, increasing the mass of the 
beamsplitter by a modest factor might already be sufficient in 
order to make the bright port fluctuations negligible for the 
GEO 600 topology. Also strategies which yield an increase of 
the sensitivity at high frequencies [10] are not influenced by 
our results. 

There is another generic mechanism existing by which a 
bright-port dark-port coupling is built up. If the transfer func- 
tion of the two arms are not equal, then the input-output rela- 
tion contains the following contribution from the bright port 



o = IO(i, n, e) + 



pr{N - E) 



\-{p 2 N + t 2 E)-W 



(29) 



where p and r denote the amplitude reflectivity and transmis- 
sivity of the beamsplitter. There are different reasons why the 
two arms are not described by the same transfer function. One 
reason could be that one arm is detuned intentionally in order 
to transmit some carrier light towards the detector where it 
serves as local oscillator for a homodyne measurement of the 
signal. For a small detuning of just one arm, the carrier light is 
transmitted into the phase quadrature of o governed by a trans- 
fer function which is proportional to the detuning. If the car- 
rier light is to be transmitted into the amplitude quadrature of 



o, then one can choose the detuning of one arm to be it/ 2 and 
the detuning of the other arm to be slightly less than n/2. In 
that case, the west port has to be detuned accordingly to main- 
tain the power-recycling condition (which cannot hold exactly 
here, since the two arms are detuned relative to each other). 
An unintentional reason for different transfer functions of the 
two arms could be that the transmissivity and the reflectivity 
of the beamsplitter are not the same. Then, radiation-pressure 
fluctuations in the two arms would be different by virtue of 
the different powers of the two respective carrier fields. Car- 
rier light is then transmitted into the phase quadrature of 5 
and the corresponding transfer function is proportional to the 
difference of the power reflectivity and transmissivity of the 
beamsplitter. There would also be a small loss of the optical 
signal due to a partial transmission into the bright port which 
is proportional to the same difference of power reflectivity and 
transmissivity. 
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